y=[3 1]
u=[8 6]

dot(u,y)
dot(u,u)
yu=dot(u,y)/dot(u,u)*u

y-yu
norm(y-yu)
